Probing exotic spin correlations by Muon Spin depolarization measurements with 

applications to spin glass dynamics. 
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We develop a new method to probe the local spin dynamic autocorrelation function, using magnetic 
field dependent muon depolarization measurements. We apply this method to /iSR experiments in 
the dilute Heisenberg spin glass AgMn(p at. %) at T > Tg where the correlations of the Mn local 
magnetic moment are strongly non-exponential. Our results clearly indicate that the dynamics of 
this spin glass cannot be described by a distribution of correlation times. Therefore, we analyze 
the data assuming a local spin correlation function which is the product of a power law times a 
cutoff function. The concentration and temperature dependence of the parameters of this function 
are determined. Our major conclusion is that in the temperature region close to Tg the correlation 
function is dominated by an algebraic relaxation term. 

PACS numbers: 



I. INTRODUCTION 

In the vicinity of a second order magnetic phase tran- 
sition, the spin-spin dynamical autocorrelation function 



(S,(<)S,(0)) = g(t), 
must take the general form 

g(t) ^ t-^fit/r) 



(1) 



(2) 



H]. In these equations i is a site index, () is a ther- 
mal and sample average, t is a time scale that limits 
the upper range of the decay, x can be defined in terms 
of the static and dynamic critical exponents, and / is a 
"cutoff" function which can take different forms. In sim- 
ple systems, such as isolated spins, the power law part 
does not exist {x — Q), and / is exponential, making q{t) 
exponential. In systems with interacting spins q{t) can 
be non-exponential. For example, in a pure ferromagnet 
close to Tc^ q{t) is given by Eq. with x > Q and ex- 
ponential /. In this case, the correlation function of all 
individual spins is the same as the global one by defini- 
tion. In complex systems, with physical properties such 
as time scale invariance [^1 or hierarchical relaxation g] 
q{t) can be strongly non-exponential. An important ex- 
ample is spin glasses with frozen-in disorder, where the 
local environments of individual spins are not identical. 
In this case numerical simulations show that the global 
correlation function also takes the form of Eq. but / is 
a "stretched exponential" [Q,||. In addition, these simu- 
lations show that Eq. E^ is valid even locally but / varies 
from spin to spin. Therefore, they have been analyzed 
with more complex cutoff functions, such as a sum of 
exponentials H , or a sum of stretched exponentials m . 



Despite this theoretical understanding of the correla- 
tion function, the complex nature of q[t) is frequently 
ignored in the analysis of experimental systems and they 
are interpreted in terms of a sum of exponential relax- 
ation rates: some spins are relaxing fast and some slowly. 
Indeed, the Laplace transform of any relaxation function 
q{t) can always be defined in terms of a spectrum of re- 
laxation times. However, the relaxation times of the dis- 
tribution are those of the modes of the system, not those 
of the individual spins. (To take a familiar analogy, in 
a pure crystal there is a spectrum of phonon modes but 
any atom is vibrating in the precisely the same complex 
way as any other atom). 

In this work we avoid the assumption of a distribu- 
tion of correlation times and determine experimentally 
q{t) using the alternative hypothesis of a unique non- 
exponential q{t). For reasons which are described in 
Sec. P our technique of choice is Muon Spin Relaxation 
[/iSR~ Muons, are implanted in a sample initially fully 
polarized, and are depolarized in the neighborhood of 
local magnetic moments. The form of the muon depolar- 
ization is a signature of q{t). The practical difficulty is to 
establish in the general case an inversion scheme that can 
take one unambiguously from the muon depolarization 
pattern back to the correlation function. Accordingly, 
our main objectives in this work are: (I) to develop (in 
Sec. Ill) a method by which one can determine the func- 



tional form of q{t) by measuring both the field and time 
dependence polarization P{H, t) of the muon; (II) to ap- 
ply this method to a real Heisenberg spin glass and to 
determine x and r (see Sec. |rv| ). 

This work is a continuation of our previous publication 
where we have demonstrated the existence of an unusual 
correlation function in the spin glass ^(;Mn(0.5 at. %) at 



T > Tg at a single temperature g. There it was shown 
that at high fields (>120 G) and late times, the muon 
polarization obeys the scaling relation 



P{H,t)^P{t/H'--), 



(3) 



and this relation was traced back to the correlation func- 
tion given by Eq. 0. However, an interpretation of the 
parameters a; or t was not provided, and no attempt was 
made to discuss the singularity of q{t) at t = 0. Here 
we extend Ref. ||] in four different directions: (I) we ad- 
dress the i = singularity, (H) we apply lower fields 
where P{H,t) is most sensitive to /, (HI) we vary the 
temperatures, and (IV) we examine three different sam- 
ples of AgMn{p at. %) with p = 0.1, 0.3, and 0.5 where 
Tg = 0.51, 1.45, and 2.80 K respectively. This allows 
us to estimate the various parameters of q{t) as a func- 
tion of p and T. Our most important finding is that q{t) 
depends on temperature mainly through x{T) in these 
samples. 
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FIG. 1. Asymmetry vs time in two samples for various lon- 
gitudinal magnetic fields. The solid lines are fit to the model 
decribed in the text. 



II. EXPERIMENTAL ASPECTS AND RAW DATA 



As demonstrated in Ref. M and derived again in 



Sec. HI a clear evidence for an exotic correlation func- 
tion is the scaling relation of Eq. ||. In order to find 
these type of scaling experimentally from the raw muon 
spectra with good certainty, one must vary the external 
field over two orders of magnitude. The /xSR technique is 
specially adapted for this task since the signal intensity 



is independent of H, and one can vary the magnetic field 
considerably, the upper limit being given by the fields 
available in the experimental set ups, and the lower limit 
being related to the range of applicability of the analysis 
method. In addition, experiment shows that the typi- 
cal time scales are of the order of a microsecond, which 
falls between the ranges covered conveniently by neutron 
techniques and by conventional susceptibility techniques. 
Our /xSR experiments were performed in ISIS where 
magnetic fields up to 4 kG are available. In these experi- 
ments we measure the time and field dependent asymme- 
try A{H, t) of the decay positrons in the direction par- 
allel and anti parallel to the initial muon polarization. 
This asymmetry is proportional to P{H, t) where H is 
applied in the direction of the initial polarization, and 
P(iJ,0) == 1; for more details see Ref. |||. The fields 
applied for each sample obey the condition fi^ffH < Tg, 
where /ie// is the magnetic moment of the Mn ions. This 
condition ensures that the impact of the fields on the lo- 
cal moment correlation function is minimal. Standard 
frequency dependent susceptibility measurements on spin 
glasses show that a magnetic field has negligible effect on 
the spin glass relaxation, even at quite low frequencies 
corresponding to times of the order of a millisecond; as 
the muSR technique is intrinsically limited to time scales 
of the order of a microsecond or less, we can safely ignore 
any direct effects of magnetic fields on the local moment 
relaxation. 
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FIG. 2. The asymmetry from Fig. la plotted as a function 
of t/H^ for fields higher than 240 G. This plot demonstrates 
the experimental validity of Eq. pi 

We have used AgMn since Ag nuclear moments are neg- 
ligible and there is no need to worry about background 
relaxation from Ag nuclei. Though dilute, the Mn nu- 
clear moments are large and it turns out that the Mn 
nuclear term can influence the muon depolarization sig- 
nal in zero applied field. The samples were prepared 
by melting in a sealed quartz tube, followed by strong 



cold work. The samples were transformed from buttons 
to sheets about 0.2 mm thick and were then annealed 
for a week under vacuum. This technique is standard 
for ensuring complete homogenization. The glass tran- 
sition temperatures were determined using a standard 
ac-susceptibility method. For the p = 0.3 and p ~ 0.1 
concentrations we were assisted by the dilution refriger- 
ator of C. Paulsen in Grenoble. 

Representative A(7J, t) measurements in the p = 0.5 
and 0.1 samples at T = 2.9 and 0.6 K respectively, and 
at various fields, are shown in Fig. [^ panels (a) and (b). 
Clearly, H affects the muon polarization very strongly 
over more than two orders of magnitude. The larger H 
the weaker the relaxation. This well known effect is called 
decoupling. The scaling relations of Eq. Q are shown to 
hold with x — > at high fields and late times for the p = 
0.5 sample in Fig. y. Our aim is to account for the rate of 
decoupling, namely, the scaling relations, quantitatively. 



III. ANALYSIS 

First we would like to emphasize the need for a new 
type of analysis for our data. In the appendix we show 
that the muon polarization is determined by the field- 
field dynamical auto correlation function 

$(i)^72^Bi(t).Bi(0)) 

(rather than by q) where B^^ is the dynamic magnetic 
field at the muon site perpendicular to the external field 
H, and 7^ — 85.16 MHz/kG is the muon gyromagnetic 
ratio. In the high fields and late time limit wc arrive at 
the relation 



where 



P{H,t)^P[t/T,iH)] 



$(t') cos{jHt')dt' 



Ti{H) J, 



(4) 



(5) 



Since in a spin glass there are no cross correlations be- 
tween spins one expects ^(t) ex q{t). Therefore, for an 
exponential correlation function q(t) with a single corre- 
lation time we write 



where 



$(t) = 2A2 exp(-i/r) 



2A2 = 7^ (Bi) 



(6) 



In this case the muon relaxation rate is given by the 
known expression 



1 



2A2r 



However, in a spin glass, the muon could stop in a 
variety of environments and can experience different in- 
stantaneous fields or correlation times. If we allow for 
a distribution of A and/or r we can obtain an average 
polarization P by 



P{H,t) = J J p{A,T)P (- 



2A^Tt 



+ {it^rHY 



dAdr 



Nevertheless, if 7^iJ ^ l/r™'", where r™™ is the short- 
est correlation time in the distribution we expect 

PiH,t)=P{t/H^) 

which is in contrast to the experimental observation de- 
picted in Fig. 0. This brings us to one of the important 
conclusions of our experiments: a distribution of correla- 
tion times could not explain our data. 

The second problem with the standard analysis 
method emerges from the "wiggles" seen in the data at 
early times and low fields. These could not be accounted 
for by Eq. Q which does not have oscillating terms. The 
source of the oscillations are those muons for which the 
transverse magnetic field has not relaxed abruptly dur- 
ing the time of one rotation. These muons will oscillate 
around the vector sum o f both inte rnal and external fields 
at a frequency u> — j^Jh^ + S^. Since at T > Tg there 

is a distribution of B"^ with zero average, the contribu- 
tion of the dynamic field will show up as a relaxation of 
the wiggles, while the frequency of the wiggles will be 
at uj = ^fiH. Roughly, if there is a cutoff time r, then 
by the time t > t abrupt field relaxation occurs for all 
muons. Therefore, wiggle pattern will not be seen for 
longer times. As a result, the observation of a "wiggle" 
pattern up to time t is the proof of local moments which 
do not relax abruptly within this time scale, and so gives 
a model independent estimate of r. The fit procedure de- 
scribed below provides a quantitative method to include 
this type of observation within the global analysis of the 
field dependent muon depolarization. 

Having demonstrated that neither the field-time scal- 
ing relations nor the wiggles in the data could be ac- 
counted for with standard analysis methods we turn to 
describe our alternative approach. If one does not take 
the long time limit that led to Eqs. andH one finds (see 
appendix) that the relation between P{H,t) and $(i) is 
via the expression 



P{H,t) = Poexp{-r{H,t)t), 



where 



r{H,t)t= / {t-t')^{t')cos{jf,Ht')dt' 
Jo 



(7) 



(8) 



TiiH) 1 + {j^tHY 



This expression can produce the oscillations seen in the 
data (for the exponential case given by Eq. see M). 



Motivated by the success of Eq. g to account for the pre- 
hminary data of Ref. ||^ we examine here the correlation 
function 



$(i) = 2A^ 



/, _/" w exp(-t/r;) 

[t + Tef 



(9) 



where Te and t; are an early and late time cutoff re- 
spectively. This ^(i) is properly normalized at i = 0. 
From neutron spin echo data |10| , t^ is of the order of 
1 picosecond, and approximately independent of temper- 
ature. Therefore, we are safe to assume that the time at 
which the muon polarization is first measured (0.1 /xsec) 
is much longer than Te. In other words, our experiment 
is in the limit of t ^ t^, and we can write 



$(i) = 2Ci(i) 



where 



and 



2_x 



C = AV, 



t(t)=i-^exp(-i/ri). 



Now we express 



r{H,t)t = 2A^T^ f {t- t')i{t') cosi-f Ht')dt'. 
Jo 



(10) 



(11) 
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FIG. 3. (a) The polarization generated using Eqs. hd, hs, 
and hj, for various frequencies lj = j^H. (b) Demonstrating 
the theoretical validity of Eq. ra for high enough values of 
id = jiiH and late times. 



To be rigorous we should allow for a distribution of 
all the parameters in the correlation function. This, of 
course, will lead us nowhere. Therefore, we assume only 
a distribution of A and so assume a unique form of q(t) 
for all local spins. This assumption appears physically 
reasonable and as we shall see it allows us to account for 
our data successfully. We note that all /iSR data to date 
were analyzed using this assumption, and, in particular, 
the following distribution |1^ was chosen 



P(A) 



1 A* 



1 



27rA2 



exp 



(12) 



As before the average muon polarization P is given by 
P{H, t)^Po f p(A) exp {-T{H, t)t) dA 



which leads to 



P{H, t) = Po exp (- [F* (H, t)t] '/' 



(13) 



where 



T*{H,t)t^2C* / {t-t')L{t')cos{jHt')dt', (14) 
Jo 



and 



C* = {A*fT^ 



(15) 



The difference between P before averaging over A, and 
P (after the average) is fundamental and is manifested in 
the appearance of a power 1/2 in the exponent of Eq. |l3|. 
In fact, the reason for choosing p{A) given by Eg . |12| 
is that experimental data agrees well with Eq. |l3|. In 
contrast, the difference between F and F* is minor; C is 
simply replaced by C* , namely, A is replaced by A*. In 
reality we do not use Eq. O but rather fit our data to 



P(i?,t) =Po exp (-[F*(i/,i)i]' 



(16) 



with /3 ~ 1/2 a global parameter for all samples and 
temperatures. 

It is instructive to examine Eq. 16 and nM numcricall 
for a particular example of t{t). The integral of Eq. |l 
is well defined and could be evaluated by the "improper 
integration" method ||l^. In Fig. ga we depict the nu- 
merical P(w, t) vs. t for various values of u; = jpH where 
the parameters of the correlation function are C* = 1, 
X = 0.3, Ti = 3.3, and in Eq. |l^ /3 = 0.5 (for other cases 
see Ref. [|3[). There are two important features in this 
figure. First, at low frequencies there are pronounced 
wiggles in the simulation waveform at early times, as in 
the data. The time scale on which wiggles are observed 



< 



Tl- 



[tw], which is demonstrated on the graph, obeys ii 
Therefore, as mentioned before, r/ could be estimated by 
visually inspecting the low field early time data. Second, 



there is absolutely no difference in the polarization be- 
tween a; = and u> = 0.3. In other words, at low fields 
the polarization is field independent and the field impacts 
the polarization only for lo > 1/t;. This provides a com- 
plementary qualitative method of estimating tj. In fact, 
both methods are well known for the pure exponential 
correlation function. The new point being made here is 
that these methods are correct even when a power law 
multiplies the exponential decay term in the expression 
for q{t), so the direct observation of the threshold value of 
the field dependence of the polarization is a complemen- 
tary, model-independent method of estimating the cutoff 
time. 



the first term we find that 
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FIG. 4. X (from Eq. hol) plotted vs. temperature for differ- 
ent samples. 

In Fig. ||b we show the numerical P{uj,t) vs. tju:^''^ 
for w > 5. Clearly, at high enough frequencies and late 
enough times scaling of the form of Eq. indeed holds. 
The critical value above which H and t are considered 
large is determined by r/. Late time is i > t;, and high 
fields are H > {jfiTi)^^. It is therefore very important 
to observe wiggles and/or strong field effect in the data 
(e.g. Fig. |l|) before the cutoff time and x can be obtained 
unambiguously from the scaling. 

We now apply these concepts to the analysis of our 
fiSK data. The analysis of the measurements shown in 
Fig. |l|(a) is discussed in detail using four steps. First, 
we estimate t;. The wiggles in the waveform prevail for 
iu) ^ 2 /isec, and fields as small as 30 G (w = 2.5 MHz) 
impact the polarization. Both these observations indicate 
that Ti is of the order of 1 /xsec. Second, we evaluate 
X. From the value of r/ we expect scaling to hold for 
t > 1 /isec and for fields bigger than 100 G. Indeed, in 
Fig. H when t/H^ > 0.01 (and even much earlier) all data 
sets collapse onto one line, meaning x — > (as mentioned 
before). Third, we estimate C* and /3. By expanding t{t) 
in powers of t, performing the integral of Eq. O term by 
term for the zero field case {H = 0), and keeping only 



r*(0,i)t 



2C*t2 



(2-x)(l-a;) 



+ 0(t''-^). 



Using this relation in Eq. |T^, and fitting the lowest field 
early time data to a stretched exponential (while holding 
X fixed from second step) we can obtain C* and /3. The 
values found up to now for C*, x, ti, and /3 are used 
as initial guesses for a global fit (at a given T) for all 
applied fields. In the global fit we allow freedom in H 
for fields < 30 G in order to account for the sample's 
susceptibility, demagnetization, and also since our power 
supply does not give accurate values in such small fields. 
The global fit produces values of x somewhat larger than 
those obtained from scaling and smaller error bars. The 
results of this fit procedure are presented in Fig. |^(a) 
by the solid lines. Our model captures the essence of 
the wiggles, although for some fields it fails to capture 
the fine details of the waveform at early times. More 
importantly, the model accounts very well for the data 
past the wiggles for more than 2.5 orders of magnitude 
in field. The same analysis method is applied to our 
series of samples at various temperatures. We find that 
/3 = 0.45(5) globally for all samples and temperatures. A 
second example of the success of this analysis procedure 
is depicted in Fig. |l|(b) by the solid lines. 

At first sight it would appear that with so many fit 
parameters no significant conclusions could be reached. 
In fact, because a single set of fit parameters are being 
used to fit a whole family of curves and because different 
parameters dominate the behavior at different ranges of 
time and field, the fit parameters are well pinned down. 
There is some play-off between different parameters so 
the values are correlated, but it turns out that the con- 
clusions that we will reach (relaxation dominated by x, 
cutoff times comparable or longer than a microsecond) 
are robust. To get fully reliable and unique fits, certain 
purely experimental parameters (background and initial 
anisotropy for the whole range of fields used) must be 
well established from independent control runs. 



IV. RESULTS 

We now turn to discuss the outcome of the fits. The 
power X is plotted vs. T/Tg in Fig. |[ At all temperatures 
studied (from Tg to about 2Tg) we find x values that are 
clearly bigger than zero, and that increase with increas- 
ing T, in a qualitative agreement with Ogielski's simula- 
tions [ ij and with the more recent results of Franzese and 
Coniglio [||. In general, at a given T/Tg, the power x ap- 
pears to increase with increasing magnetic ion concentra- 
tion p. Theoretically, near a continuous phase transition 
the power x is expected to relate to the static and dy- 
namic critical exponents through x = {d + jf)/2z H]. We 
find X values which are of the order of 0.15 close to Tg. 
This should be compared with the value of 0.13 that can 



be estimated from exponents obtained by magnetic mea- 
surements at low frequencies on the same type of alloy 
Q . This exponent could be expected to be independent 
of concentration for dilute alloys. 

As for the cutoff time r/, we have noticed that r/ 
changes little as a function of the temperature over the 
temperature range we have studied. This is best demon- 
strated in the p = 0.1 sample where the waveform has 
strongest wiggles and strong field dependence, and is 
therefore most sensitive to rj. In Fig. ^(a) and (b) we 
present P{H,t) at three representative fields and two 
different temperatures approaching Tg. Clearly a field 
as small as 5 G (w = 0.4 MHz) impacts the polarization 
at both temperatures, and there is no big change in the 
wiggle time tu, (~ 10 /isec) between these temperatures. 
Thus Ti remains longer than a few microseconds up to 
T/Tg = 1.43. In the inset of Fig. |(a) we show the fit 
results for 1/r; for the p = 0.1 sample, demonstrating 
once again that we do not observe a critical temperature 
dependence. On general physical principles related to 
any continuous transition, t; should diverge as T tends 
to the ordering temperature (see for instance Ogielski's 
Ising spin glass simulations [H). The present data sim- 
ply indicate that for the concentrations and temperature 
ranges studied, n is at the edge or out of the muon dy- 
namic window. 
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FIG. 5. Demonstrating that neither the wiggling time t^, 
nor the sensitivity of the polarization to the field are varying 
with temperature. The inset shows the late cutoff time r; as 
given by Eq. hol versus T/Tg for the p = 0.1 sample. 



be T independent, and so should (A*)^. If Tg is also T 
independent, we expect from Eq. |l^ a log(C*) against x 
plot to be a straight line. Indeed, a linear dependency be- 
tween log(C*) and x is consistent with the data as shown 
in Fig. 0. The slope of the line at each concentration is 
log Te and the intercept is 2 log(A*). The values of the in- 
tercepts have large error bars (since we do not have data 
in sufhciently small x) and are therefore not presented. 
Nevertheless, a rough estimate of A* indicates that the 
condition for validity of our analysis, namely, H > A* 
(see appendix) holds for H ^ 10 G for the p = 0.5 and 
0.3 samples, and H ^ 1 G for the p > 0.1 sample as 
done here. The values of logTg, are presented in the in- 
set of Fig. ^ for the different samples. These values make 
the conjecture t ^ Te self-consistent, although we find Te 
for the p = 0.1 sample unacceptably small. Nevertheless, 
our analysis indicates that the T dependence of C* is due 
only to x{T). 



t_< 



lE-3 
O 




O.1°/o 



. 1 0.2 0.3 0.4 <).5 O.f) 0.7 0.8 



FIG. 6. The prefactor C* plotted versus a; on a semi-log 
scale (see Eq. |l5[ ). The temperature is an implicit parameter. 
The solid curves are linear fits. Their slope is interpretted as 
logTe and presented in the inset. 

To summarize, we demonstrate that a distribution of 
correlation times cannot account for the muon polariza- 
tion as a function of both field and time. However, we 
can account for our data provided that each muon expe- 
rience a local field correlation function given by Eq. |9|, 
where the only difference between different muons is in 
the value of A. Our analysis shows that out of all the 
parameters that determine P{H, t), namely, x, ti, Te, and 
A*, the only one that is temperature dependent is x. We 
thus conclude that the dynamic behavior of this family 
of dilute spin glasses near Tg is controlled essentially by 
the temperature dependence of the power law exponent 

X. 



Finally, we examine the parameter C*. The mean 
square of the field (B^) (or A^) at a given site should 
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APPENDIX A: FROM FIELD CORRELATIONS 
TO SPIN LATTICE RELAXATION 

In this appendix we derive the expected behavior of 
the polarization of a local probe, which is out of equilib- 
rium, in a dynamic field environment. For simplicity we 
use a semi-classical approach. A more complete result, 
based on a full quantum treatment, can be found in the 
literature [15[ but there is no detectable difference. In the 
semi-classical method we decompose the Hamiltonian of 
a spin 1/2 experiencing both an external static field H, 
and an internal dynamically fluctuating field B'*(t) into 

n = na + n'{t) (ai) 

where 

Wo = -J^.S,H (A2) 

is the secular part, and the interaction part is 

H' = -7mS • B^{t). (A3) 

The time dependent field B'^(i) is taken to be classical 
and S is the muon spin operator. When the fluctuat- 
ing flelds are smaller than the external fleld we can use 
time dependent perturbation theory and write the time 
propagator as 

Uit) = exp(- jHot)[l -^J dt'n\t') - (A4) 

^ j dt' I dt"n\t')n\t") + ..]. 

where the perturbation Hamiltonian in the interaction 
picture is given by 

H^{t) = cxp{mot/h)n'{t) e^i>{~inot/h). (as) 

This Hamiltonian simplifles to 

n\t) = ~j,B\t)S\t) (A6) 

where 

S (i) = exp{iujSzt/Ti)Sexp{—iujSzt/h) 
and 



SKt) = SxCOs{ujt) — SySin{Ljt) 
SUt) = SyCos{ujt) + S^sm{Ljt) 

slit) = Sz. 



(A7) 



Note that S^(i) is the time-dependent spin operator in 
the interaction picture, namely, of a muon that rotates 
around the external fleld as if there were no internal 
fields. It is not the spin operator S(i) of a muon that 
experiences the combined static and dynamic fields. In 
other words, in order to find S(f) we use a perturbation 
expansion in terms of S^(i). This is the (mathematical) 
origin of the wiggles seen in the data at the frequency of 
the external field. 



Using Eqs. A6 and A7 we find that 

n\t) = -7^ [V{t)S^ - U(t)Sy + T{t)S,] (A8) 

where 

Vit) = B^{t) cos{Lut) + B^{t) sin{ujt) 
U{t) = B{{t)sin{ujt)-Bl(t)coslujt) 
T{t) = B^^it). 

The polarization of a muon at a given site as a function 
of time Pz (i) is given by 



Pzit)^TT[pU\t)azU{t)] 



(A9) 



where 



1 + Po<yz 



and Pq is the initial polarization. 

Equations A£ and A4 lead to the perturbation series 



P{t) =1 + 11 + 111 



(AlO) 



where 



/ = Ti- {pa z}, 

II = fJ*dt'TTp[H'{t'),az], 

III = -^ /o dt' J^' dt"Trp[H\n, [n'{t'),az]]. 

(All) 

We now evaluate each term explicitly: The first term is 
simple. Since 



TrIpcrJ = 



Pn 



i = x,y 

we get I = Pq. In the second term we have 

[n\t'), C7z] = z7m iVit'Wy + U{t')a,) . 

After the evaluation of the trace we find II = 0. For the 
third term we find 

Trp[W'(t"), [H'{t'),az]] = Poll [V{nV{t') + C/(t")C/(i')] 



where we have used the simplifying assumption that V, 
f/, and T are classical fields. To further progress with 
the calculation it is helpful to introduce at this stage the 
ensemble average over all possible trajectories of B'[(i). 
Quantities in () will denote such an average and we can 
write 



where (/) — Pq and 

(///) ^ Po^l {V{nV{t') + U{t")U{t')) . 



(A12) 



The meaning {Bf{t')B'^(t")) is to hold t' and t" fixed 
and to average over all possible fields in the systems in 
between these times. This is a correlation function. Wc 
assume that there are no cross correlations, and x,y sym- 
metry so the only terms left to evaluate are of the form 

{V{t")V{t') + U{t")U{t')) 
= {Bi{t')Bi{t") + B^y{t')B^y{t")) cos {u{t" - t')) . 

We now make the assumption that the correlation func- 
tion depends only on the time difference and define 

W - t") = ll {Bi{t')Bi{t") -f B^^{t')B^^{t")) . (A13) 

so that finally 

(///) = -Po / dt' f dt"^{t' - t") cos {uj{t' ~ t")) . 
Jo Jo 

We can now eliminate one of the integrals by writing 

t" 



T ^ t' 
t' ^ t- 



t' 



This transforms the integral into 



dt' I dt" 

^0 



dr 



dr' 



and we arrive at 



(///} = -Po f dT{t - t)$ 

Jo 



(r) cos(a;T). 



If we now write 

P^(i) = Poexp(-r(i)i), 



(A14) 



(A15) 



expa nd th is e quatio n in powers of F, and compare it with 
Eqs. |A12| and A14| we find 



T{t)t= I dT{t-T)<P{T)cos{ujT). 
Jo 

At late times, such that <i>(i) is negligible, one finds that 
1 



r(ij,i) 



Ti{H) 



^(t) cos(cjr)dr. (A16) 



[1] P. C. Hohenberg and B. L. Halperin, Rev. Mod. Phys. 

49, 435 (1977). 
[2] H. Scher, M. F. Shlesinger, and J. Bendler, Physics To- 
day, January, 26 (1991); J. Klafter, M. F. Shlesinger, and 
G. Zumofen, Physics Today, February, 33 (1996) 
[3] R. G. Palmer, D. L. Stein, E. Abrahams, and P. W. An- 
derson, Phys. Rev. Lett 53, 958 (1984). 
[4] T. Ogielski, Phys. Rev. B 32, 7384 (1985). 
[5] G. Franzese and A. Coniglio, Phys. Rev. E 59, 6409 

(1999). 
[6] H. Takayama and H. Yoshino, Physica A 204, 650 (1994). 
Hiroshi Takano and Seiji Miyashita, J. Phys. Soc. Jpn. 
64, 560 (1995). 
[7] S. C. Glotzer, N. Jan, T. Lookman, A. B. Maclsaac, and 

P. H. Poole, Phys. Rev. E 57, 7350 (1998). 
[8] A. Keren, P. Mendels, L A. Campbell, and J. Lord, Phys. 

Rev. Lett. 77, 1386 (1996). 
[9] A. Keren, Phys. Rev. B 50, 10039 (1994). 
[10] F. Mezei and A. P. Murani, JMMM 14, 211 (1979). 
[11] Y. J. Uemura et al, Phys. Rev. B 31, 546 (1985). 
[12] W. H. Press, B. P. Flannery, A. A. Teukolsky, and W. 
T. Vetterling, Numerical Recipes (Cambridge University 
Press, Cambridge, 1989). 
[13] A. Keren and G. Bazalitsky Physica B 289, 205 (2000). 
[14] L. P. Levy and A. T. Ogielski, Phys. Rev. Lett 57, 3288 

(1986); L. P. Levy, Phys. Rev. B 38, 4963 (1988). 
[15] T. McMulIen and E. Zaremba, Phys. Rev. B 18, 3026 
(1978). 



